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Abstract 

Partial pair-correlation functions of colloidal suspensions with continuous polydispersity can 
be challenging to characterize from optical microscopy or computer simulation data due to in- 
adequate sampling. As a result, it is common to adopt an effective one-component description 
of the structure that ignores the differences between particle types. Unfortunately, whether this 
kind of simplified description preserves or averages out information important for understand- 
ing the behavior of the fluid depends on the degree of polydispersity and can be difficult to 
assess, especially when the corresponding multicomponent description of the pair correlations 
is unavailable for comparison. Here, we present a computer simulation study that examines 
the implications of adopting an effective one-component structural description of a polydisperse 
fluid. The square-well model that we investigate mimics key aspects of the experimental behavior 
of suspended colloids with short-range, polymer-mediated attractions. To characterize the par- 
tial pair-correlation functions and thermodynamic excess entropy of this system, we introduce a 
Monte Carlo sampling strategy appropriate for fluids with a large number of pseudo-components. 
The data from our simulations at high particle concentrations, as well as exact theoretical results 
for dilute systems, show how qualitatively different trends between structural order and particle 
attractions emerge from the multicomponent and effective one-component treatments, even with 
systems characterized by moderate polydispersity. We examine consequences of these differences 
for excess-entropy based scalings of shear viscosity, and we discuss how use of the multicom- 
ponent treatment reveals similarities between the corresponding dynamic scaling behaviors of 
attractive colloids and liquid water that the effective one-component analysis does not capture. 



I. INTRODUCTION 

Fluid suspensions of natural and technological interest are inevitably polydisperse. 
Their constituent particles-each of which comprise a large number of smaller atoms, molec- 
ular, or ions-exhibit distinct geometric and chemical features (size, shape, charge, etc.) 
and hence represent different species of a many-component mixture. As might be expected, 
the composition of such a mixture can strongly affect its structural, thermodynamic, and 
dynamic behavior.-"— For example, suspensions with narrowly distributed (e.g., approx- 
imately monodisperse) particle characteristics often readily crystallize and thus can be 
interesting candidates for applications like colloidal self assembly.— "^I In contrast, fluids 
with more diverse particle populations that suppress crystallization:^»ii»i^'^i^i^»2^'^i2^ are 
more useful in applications which call for good glass-forming materials. 

Characterizing the properties of systems with a large number of components can be a 
formidable technical challenge. Insufficient sampling typically limits the species-specific 
information that can be directly obtained via experiments and computer simulations. In 
fact, for suspensions with apparently "weak" polydispersity, it is common practice to ig- 
nore particle differences altogether, adopting an effective one-component description. The 
latter is tantamount to assuming a fictitious pure-fiuid model, whose (typically softer) 
interparticle interactions are those that, under equilibrium conditions, would reproduce 
some globally averaged behavior of the original polydisperse material.— i^"^"— An obvious 
drawback of this simplification is that one generally does not know in advance whether a 
system's polydispersity is weak enough (as it pertains to the properties of interest) to jus- 
tify use of an average description that washes out species-specific contributions. Moreover, 
it is difficult to assess the impact of averaging a posteriori, especially when statistically 
meaningful information about the component contributions to the polydisperse material's 
properties is unavailable. This type of dilemma, though far from resolved, has long been 
appreciated and studied in colloid science. For example, it arises in considering the trade- 
offs between characterizing polydisperse colloidal suspensions by partial (species-specific) 
versus total (average) structure factors inferred from scattering data.-'^'^'^'^'^'^""— Al- 
though a multicomponent description is ideally desired, it cannot generally be obtained 
from experiments without invoking a number of approximations that compromise its reli- 
ability. 



In this study, we use computer simulations and exact theoretical results to examine 
how the static structure of a model polydisperse fluid differs when quantified based on 
multicomponent versus effective one-component descriptions of its pair correlations. Sim- 
ilar in spirit to, e.g., related work of Pangburn and Bevan,— i^ we focus exclusively on 
the real-space counterpart to the reciprocal-space scattering problem discussed above, the 
former being especially critical for the interpretation of data from computer simulations 
and optical microscopy experiments^"— of complex fluids. 

The model we investigate here is a fluid of polydisperse hard-sphere particles with 
short-range, square-well attractions. Earlier simulations^ show that it reproduces key 
static and dynamic behavior of suspended colloids with polymer-mediated depletion in- 
teractions. Most notably, in agreement with experiments^'^ it displays strikingly anoma- 
lous dynamic behaviors at high particle concentrations near the so-called "repulsive" glass 
transition (e.g., diffusivity increases and viscosity decreases with increasing particle at- 
tractions). Interestingly, an effective one-component analysis^ suggests these behaviors 
emerge under conditions where the fluid also displays structural anomalies (e.g., static 
correlations apparently weaken with increasing particle attractions). Of course, it is im- 
portant to note that polydispersity is an important feature of this and similar models 
because it suppresses crystallization that would otherwise restrict access to the deeply 
supercooled fluid states where the anomalous behavior is most pronounced. 

Here, we use exact results to probe some of the consequences of adopting the effective 
one-component description for the structure of the aforementioned model in the dilute 
limit. To carry out a similar analysis at high particle concentrations, we introduce Monte 
Carlo sampling techniques that allow us to accurately determine the m{m + l)/2 partial 
radial distribution functions as well as the free energy for a m pseudo-component [in this 
case, m = 60; m{m + l)/2 = 1830] description of the model. 

The main findings of this study are as follows. (1) The multicomponent (i.e., partial) 
and effective one-component radial distribution functions show qualitatively different be- 
havior as a function of the magnitude of the interparticle attractions, even for this moder- 
ately polydisperse model. This serves as a cautionary note regarding conclusions that one 
might otherwise draw about structure in polydisperse systems based on trends exhibited 
by the latter globally-averaged, static correlations. (2) The pair-correlation contribution 



to the excess entropy of the polydisperse fluid based on the multicomponent description, 
in contrast to the estimate based on the one-component model, very closely approximates 
the full thermodynamic excess entropy of the system. This suggests that techniques which 
allow one to accurately determine the partial radial distribution functions of dense, poly- 
disperse colloids will also furnish information key for understanding their thermodynamic 
behavior. (3) Finally, the excess entropy scaling for the shear viscosity of this polydis- 
perse fluid is similar in some key ways to that of liquid water, another system known to 
exhibit anomalous thermodynamic and dynamic properties. As we show, such similarities 
are not apparent if the effective one-component description is used in the excess entropy 
estimation for the polydisperse colloidal fluid. 

II. METHODS 

A. Model for the polydisperse fluid 

We consider a previously introduced model^- comprising hard-sphere particles with 
short-range, square-well (SW) attractions. The pair potential between particles i and j 
(with diameters Uj and aj, respectively) is given by 
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Here, r is the separation between particle centers, ctjj = (uj + crj)/2 is the hard-core con- 
tact separation, and e is the depth of the attractive well. The particle diameters o"j of 
this polydisperse fluid follow a truncated normal distribution with a mean value of a and 
a standard deviation of O.la. The truncation ensures that all particle diameters of the 
fluid are within 3 standard deviations of the mean (0.7cr — 1.3o"). This degree of polydis- 
persity in the particle population, while still modest, is signiflcant enough to discourage 
crystallization in simulations, allowing access to deeply supercooled fluid states.— 

As discussed extensively elsewhere,— i^i^ short-range attractions like those exhibited by 
this model can give rise to unusual dependencies of shear viscosity (which decreases) and 
self diffusivity (which increases) with e/k^T, where /cb is the Boltzmann constant and T 



is temperature. These trends capture what is observed experimentally^'^ for suspensions 
of hard-sphere-like colloids with polymer-mediated depletion attractions. Interestingly, an 
effective one-component analysis of this model's pair correlations^ suggests that the afore- 
mentioned dynamic anomalies are accompanied by a structural anomaly (i.e., increasing 
e/ksT weakens the effective pair correlations). In this work, we test whether this feature 
also emerges from a multicomponent analysis. 

B. Monte Carlo simulations 

For the Monte Carlo (MC) simulations, we approximate the model introduced in Sec- 
tion III Al by a closely related mixture of 60 pseudo-components. Particle diameters in the 
pseudo-component mixture exhibit continuous polydispersity, and their pairwise interac- 
tions Vij{r) are described by eq. [1] The pseudo-components each comprise nearly identical 
particles (e.g., diameter differences of less than O.Olo"). Specifically, particles with diame- 
ters greater than [0. 7-1-0. 01(A; — l)]cr and less than [0.7-1- O.OlfcJa are considered members of 
the k^^ pseudo-component with chemical potential /i^, where k G [1,60]. In what follows, 
subscripts on variables refer to pseudo-component numbers rather than particle identities. 
We have found that specifying the nominal mole fraction of the k^^ pseudo-component, 
x^, by the relation 

„ ^ exp[- (A: -30.5)7 200] ..^ 

""' Ef=iexp[-(/-30.5)V200] ^^ 

results in a mixture with properties virtually identical to the model described in Sec- 
tion III A[ As we discuss below, the pseudo-component mixture is very convenient to 
simulate within the semi- grand (SG) ensemble.— 

In our SG MC simulations of this mixture, temperature T, total particle number A^ = 
1000, volume V (or density p = N/V), and all pseudo-component chemical potential 
differences Aprk = Pk — Pr (relative to that of an arbitrarily chosen r*^ pseudo-component 
reference) are held constant. In what follows, we express chemical potential differences 
in terms of activity ratios ^rk = C,k/C,r, where the activity of a particle of the fc*^ pseudo- 
component with molecular partition function g^ is given by ^k = qke^''^''^^. 

In order to study a fluid with a specified pseudo- component composition x^ within the 



SG ensemble, we must first solve for the activity distribution i^°j. that produces the desired 
x'^ at the temperature and density of interest. In general, this is a highly nontrivial task 
because ^°^ is an unknown functional of x^. Fortunately, this inverse problem has received 
considerable attention,—"— and robust methods are now available for computing Qf^ given 
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Here, we use an efficient nonequilibrium potential refinement scheme introduced by 
Wilding^ to obtain ^°^. We take advantage of the efficient sampling of composition space 
facilitated by SG MG simulation techniques and allow particles to adopt any diame- 
ter within the range 0.7a — l.Scx. We begin the potential refinement scheme by setting 
In^rfc = and subsequently adjust the activity distribution during an eight-stage SG MG 
simulation. At regular intervals during the i*^ stage, In^^^ is incremented by the relative 
difference between the instantaneous Xk and target x^ discretized particle-size distributions 
scaled by a modification factor 7j, 

lnCr = lnCf-7.(^^). (3) 

A given stage terminates when the maximum relative difference ( = max [|(a;^^^ — x^) /x^|] 
between the target (x^) and stage-averaged (x^^^) particle size distributions drops below 
a tolerance of (* = 0.01. We set 71 = 0.001 and reduce the modification factor by a factor 
of two after the completion of each stage. The activity distribution that emerges after the 
eighth stage is taken to be ^°^. 

C. Static pair correlations 

Perhaps the most basic nontrivial measure of real-space structure in a homogeneous, 
multi-component mixture is the partial radial distribution function (PRDF), gij{r). The 
PRDF characterizes the frequency of various pair separations r that occur between parti- 
cles of "type" i and j in the fiuid. Specifically, the mean number of particle centers of type 
i located in a spherical shell of differential thickness dr a distance r away from a type j 
particle center is 4TTr'^pigij{r)dr, where pi is the overall number density of particle type i. 
Note that symmetry requires gij{r) = Qjiir). Below, we use the labels i and j to indicate 
specific choices of the m = 60 pseudo-components, i.e., z,j G [1,60]. Note that there are 
m{m -|- l)/2 = 1830 distinct pseudo-component PRDFs in the mixture we study here. 



highlighting the statistical challenge associated with characterizing the structure, even at 
the pseudo-component level of description. 

The effective one-component treatment for pair correlations, on the other hand, ignores 
all differences between the various particle types. Its radial distribution function gcs{i^) 
is defined such that the mean total number of particle centers located in a spherical shell 
of differential thickness dr a distance r away from other centers is 47ir'^pges{'r)dr, where p 
is the total particle number density. Expressed differently, ges{r) = J2i=i J2j=i^i^j9iji^)- 
In other words, as should be expected, the effective one-component description provides 
a highly averaged representation of the fluid's pair correlations. This makes it easier to 
characterize statistically, but more challenging to interpret. 

D. Structural order metrics 

We examine the behavior of two different structural order metrics that characterize the 
degree of local translational order of a fluid based on the strength of its pair correlations. 
To understand the implications of adopting the effective one-component description for 
our polydisperse model, we compare the value that these scalar parameters take on when 
the structure is described by ges{r) versus the full set of pseudo-component PRDFs, i.e., 
gijir), where i,j e [1,60]. 

The first structural metric we consider is —s^^'/kBr^ where s^'^' comes from the leading 
term in the multiparticle expansion of the fluid's molar excess entropy (over ideal gas), 
gex _ g{2) _|_ g(3) _|_ . . . .55i56 rpj^g explicit Connection between — s'-^^^b and the PRDFs of a 
fluid mixture can be seen when the former is expressed as 

- y^ = E E ^ /t^^^^^) ^^^a^Ar) - ^,,(r) + l]dr (4) 

i J 

Note that — s'-^-'/^b is a non- negative quantity that vanishes for an ideal gas and is consid- 
erably larger for dense liquids and glasses that exhibit stronger interparticle correlations.— 
Since 3- and higher-body correlations are challenging to characterize in a fluid, the thermo- 
dynamic quantity s^^ is commonly approximated by s^^-* for pure fluids or binary mixtures. 
Here, we compute s'^^ and s^^-* (at the pseudo-component level) directly from simulation. 
The corresponding structural order metric in the effective one-component description is 
given by 
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(2) p 

- ^ = f y [^eff(r) In^efr(r) - ^eff(r) + l]dr (5) 

As discussed earlier, it was this quantity that was used in previous computer simulation 

studies to characterize the average strength of the interparticle correlations for the system 

considered here^ and related models of polydisperse colloids^. 

The second structural order metric r studied here is defined as 

rvc 
-r = Vc^ y] y] Xi^j / \9^j{y) - ^dy (6) 

where y = rp^l'^ and y^ is a cut-off value (here, we chose y^ = 4). This measure is a 

straightforward generalization of a parameter introduced earlier to study the local trans- 

lational order of fluid, glassy, and crystalline states of single-component materials based on 

their static pair correlations.— i^ The version of this metric for the effective one-component 

model is given by 

rvc 

TcflF = Vc^ / \ges{y) - Mdy (7) 

where, as before, y = rp^l^ and y^ = 4. 

E. Thermodynamic excess entropy 

A two-step process is used to compute the excess entropy per particle s^"" (pcr^ , e/k^T) 
of our pseudo-component SW fluid. In the first step, we obtain thermodynamic properties 
as a function of po"^ for polydisperse hard-spheres {e/k-BT = 0) in contact with a reservoir 
of particles with an ideal-gas activity distribution ^^| = a;^. In the second step, we move 
along a path that takes the system from ^^| and e/k^T = to Q^^ and nonzero e/k^T at 
constant density. 

To obtain the density dependence of the entropy, we perform a multicomponent grand 
canonical (GC) simulation with V = lOOOcr^ and C,k = ^rCA^ with ^r = Qr- Following the 
strategy outlined above, we allow the individual particle diameters to take on values that 
span the range 0.7a to 1.3a and treat C,k as a stepwise function. Transition matrix MC 
methods^"^ are used to determine the probability nGc(A^) of observing the system with 
a total number of particles N. The density probability distribution is related to the SG 



partition function T and GC partition function S as follows^ 

^[ij,V,T) \qrj 
The relevant bridge equation Y = —k^T InT provides the SG potential 



Y (iV) = -kj,T 



1 ,noc(iV)V^^y6. 



(9) 



Hgc (0) ; \q^ 

where we have used the zero-particle limit to express the GC partition function in terms 
of the particle number probability distribution, InS = —In 11(0). The SG potential and 
Helmholtz free energy F are linked by the relationship^ 

•m 

F{N) = Y{N) + nJ2 ^k {N) {fik - f^r) , (10) 

fc=l 

where Xk (N) is the ensemble-averaged discretized particle size distribution within a system 
described by T, V, N, and ^rk- For the e/k^T = case, the intensive entropy s and 
Helmholtz free energy / are now 



kB kuT N \Ugc{0) 






Within the GC simulations employed here , we have ^^./^ = x^, and therefore Xk (N) 
deviates from the normal particle size distribution. For the polydisperse hard-sphere fluid 
{e/kBT = 0), increasing A^ causes the Xk (N) distribution to shift such that the populations 
of pseudo-components with smaller particle diameters increase. This feature of the system 
allows us to sample high particle number densities with relative ease. To summarize, 
within the first stage of the two-step scheme employed here, we use GC simulation to 
obtain the density dependence of thermodynamic properties for a hard-sphere polydisperse 
fluid {e/k^T = 0) subject to ^^./^ = x^, and we perform just one of these simulations. 

In the second step of our approach an expanded ensemble (EE) MC procedure^ is used 
to determine the difference in the thermodynamic properties of a polydisperse SW fluid 
with C,rf^ and the reference polydisperse hard-sphere fluid with ^^^ at a pa^ and e/k^T of 
interest. The two systems are connected through a series of subensembles described by 
A, which spans from A = to 1 in increments of AA = 0.001 (i.e., Aj = O.OOli, where 
i G [0, 1000]). To obtain the relative entropy of the system along this path, we perform SG 
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simulations with fixed N = 1000 and pa^ = Na^/V. Within subensemble i, the activity 
distribution C,rk,i and well depth ei/k-sT are 

Ine.., = InCl. + A. (inC, - InCl) (12) 

and 

ei/k^T = Kt/k^T. (13) 

Transition matrix MC methods^i^ are used to evaluate the probability Hee (A) of finding 
the system in a subensemble defined by A. The difference in the SG potential over the 
entire path evaluates to 

AF = r (C\„ e/fceT) - Y (Ct O) = -k^T In (^^||^^) • (14) 

The corresponding change in the Helmholtz free energy per particle is 



A/(po-3,eABT) 



A^ 



K2si)--|:^^mKf)-^(o)'" 



tig 

Qk 



where Xk (0) and Xk (1) are the ensemble-averaged discretized particle size distributions 
within the A = and 1 subensembles, respectively. By further considering the ensemble- 
averaged intensive configurational energies e(A), we arrive at the entropy difference As 

^i^^^!^^i/^ = ^Hl)-e-(0)-A/(p.»..AeT)]. (16) 

Kb k-qI 

where, in our case e (0) = 0. An EE simulation is completed for each pa^ and e/k^T of 
interest. Collectively, the GC and EE simulations provide the absolute entropy 

s {pa\ e/k^T) = s {pa\ O) + As (pa^ e/fceT) . (17) 

We now turn our attention to the ideal gas contribution to the excess entropy. The 
entropy of a multicomponent ideal gas is 



m ^ X 

1 - In p - V" Xfc In ( — ) 



5lg _ . 



To be consistent with the analysis above, we evaluate the expression above with Xk = 
Xk (1), which, of course, closely approximates x^. Finally, we arrive at the excess entropy 

s*^- = s - s'^. (19) 

We note that the algorithm outlined here is rather general and we expect it to prove useful 
for computing the free energy of a wide range of polydisperse fluids. 

11 



III. RESULTS AND DISCUSSION 

The first results that we consider pertain to the behavior of a dilute (p — )■ 0) polydis- 
perse fluid of SW particles that interact through the potential given by eq. [1] and that 
exhibit a Gaussian distribution of particle diameters (with mean a and standard devia- 
tion cr/10). This is an interesting case to study because-while simple enough to allow for 
an exact theoretical description of the polydisperse fluid's static properties-it is also rich 
enough to illustrate some of the key consequences of adopting the effective one-component 
description. Analytical expressions for the structural quantities discussed below can be 
found in the Appendix. 

Fig. [1] shows how the strength of the static interparticle correlations of the dilute fluid 
depends on the magnitude of the interparticle attraction e/k-QT. Two levels of structural 
characterization are considered: one that takes into account the full set of PRDFs (i.e., the 
exact multicomponent description), and one based on the effective one-component struc- 
ture. The key point of Fig. [1^ is that the structural order metrics that properly account for 
the polydisperse nature of the fluid {—s^'^^/pa^k^ and r/ p^^'^a) show qualitatively differ- 
ent behavior than their effective one-component counterparts {—s\^ / pa^k-Q and Tcs/ p^ cr). 
The former display the expected monotonic behavior; i.e., static correlations strengthen 
with increasing magnitude of the interparticle attractions. In fact, in the Appendix, we 
demonstrate that this monotonic trend must always be observed in the dilute limit for this 
system, independent of the width of the diameter distribution. In contrast, the effective 
one-component quantities both spuriously indicate the presence of a structural anomaly; 
i.e., conditions where static correlations apparently weaken with increasing interparticle 
attractions. 

Why do — Sgg/pa^/cB and Tefj/p^/^a predict qualitatively different behavior than 
—s^'^^/pa^k^ and t/ p^^'^a, respectively? The basic mathematical reason is that the struc- 
tural metrics are nonlinear functionals of the PRDFs. Thus, their behaviors can depend 
sensitively on whether the latter "inputs" are treated in a multicomponent or a pre- 
averaged (effective one-component) manner. From a physical perspective, as illustrated in 
Fig. [lb, averaging the PRDFs to obtain geE{f) softens the apparent pair correlations in the 
fluid, as well as the effective interactions between the particles, Veff(r) = — kBT In gcs{r). 
As a result, the effective one-component metrics significantly underestimate the effect that 
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the attractions have on static structure. For example, the magnitude of — Sg^/pcr^/cB is 
much smaller than that of —s^'^ypa^k^ (which is equal to —s^^/pa^k^ in this dilute limit) 
for e/k-QT > 1, illustrating that one cannot readily estimate the thermodynamic excess 
entropy from the effective one-component approximation. 

In order to compare this behavior to that of concentrated SW fluids, we analyze results 
from our MC simulations. Fig. |2] displays 60 of the 1830 PRDFs that we collected for the 
polydisperse SW fluid described in Section |Tl] at a reduced density pa^ = 1.05 (particle 
packing fraction of 0.567) and two extreme levels of interparticle attractions: e/k^T = 
(top panel of Fig. |2]) and t/k-QT = 2.5 (bottom panel of Fig. [2]). The functions presented 
in Fig. |2l labeled g30j{r), characterize the partial pair correlations between particles with 
diameters close to the median size (the 30*'' pseudo-component) and particles of other 
sizes (the j*'* pseudo-component with j G [1,60]). In the case of e/k^T = (i.e., the 
hard-sphere limit), there is a clear relationship between the contact value g?Mj{r = ctsoj) 
and the contact diameter a^oj = (o"3o + <^j)/'^'i larger particles have higher probabilities 
of contact configurations. This is the expected trend and is a classic manifestation of 
the depletion effect.—'— Contact configurations of big spheres are entropically favorable 
because they increase the volume accessible to smaller spheres. A weak splitting of the 
second peak of the PRDFs is also apparent-as is commonly observed for dense hard- 
sphere fiuids^^the precise shape of which depends on the identities of the participating 
pseudo-components. Interestingly, for e/k^T = 2.5, the aforementioned entropic packing 
preferences have all but vanished. The strong attractions in this limit dominate the pair 
correlations, and there are only very minor differences between PRDFs involving different 
pseudo-components. The main point is that the sensitivity of the distribution of PRDFs to 
e/k^T suggests that the pre-averaging required for the effective one-component description 
may also be problematic for this high-density system, as it was for the dilute case. 

The data of Fig. [3] clarify this point. Specifically, panels [3K-|3]F illustrate how 5 charac- 
teristic PRDFs for this system-(730j(r) for j = 10, 20, 30, 40, and 50-evolve with increasing 
strength of the interparticle attraction. The pronounced changes of these functions [e.g., 
the contact values gmjij = ctsoj) triple in magnitude with increasing e/k-oT] should be 
contrasted with the behavior of the corresponding effective one-component quantity ges{f) 
shown in Figure |3^, which not only looks qualitatively different (considerably softer than 
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the PRDFs) but also is remarkably insensitive to the value of e/k-^T considered. In other 
words, the averaging involved in forming ges{f^) has the effect of hiding some significant 
changes that attractions have on the static structure of the polydisperse fluid. Similar-if 
less dramatic-effects emerge when adopting coarse-grained structural descriptions of other 
types of materials (e.g., center-of-mass versus monomer-level pair-correlations in fluids of 
chain-like molecules^). 

The consequences of the aforementioned averaging become clearer in Fig. H] where we 
directly compare how the structural order metrics-based on multi-component versus one- 
component descriptions of the systems-depend on e/k^T. The first point to note is that, 
as in the dilute fluid, the metrics —s^'^^/k^ and —sl^/k^ show different behavior, both 
quantitatively and qualitatively. The former provides an excellent approximation to the 
thermodynamic excess entropy —s^^/k-Q, while the latter is too small by a factor of be- 
tween 1.8 and 3.5 (depending on e/k^T). Also, as observed in the dilute case, the effective 
one-component quantity reports that the fluid is structurally anomalous (attractions ap- 
parently weaken pair correlations) across a wide range of e/k^T. In contrast, —s^'^^/k^ 
shows normal behavior (attractions strengthen correlations) except for a very narrow range 
of conditions at low values of e/k-^T where a weak structural anomaly is present. The 
quantity —s'^^/ks displays normal behavior for all e/k-QT . Similar qualitative trends are 
apparent when comparing r and Teff. The latter underestimates the magnitude of the 
structural order, and it overestimates the range of conditions where structural anomalies 
occur. 

As a final point, we examine the implications of the one-component analysis for scaling 
relationships between static structure and dynamics. In particular, we plot in Fig. [5] re- 
duced shear viscosity r]p"'^^^(mkBT)^^^'^ (here i] shear viscosity and m is particle mass) for 
the polydisperse SW fluid as a function of the various structural order metrics analyzed 
above. The shear viscosity data we use for this plot was extracted from the earlier molecu- 
lar dynamics investigation of Krekelberg et al.— . As Rosenfeld^i^ originally observed, the 
non-dimensionalized form 'r]p^^'^(mkBT)^^^^ is strictly a single- valued function of — s'^^/Zcb 
for a fluid of classical particles interacting via a pair potential of the inverse-power-law 
form and obeying Newton's equations of motion. For fluids with more complex structures 
and interactions, the scaling can only be expected to approximately hold. 
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The key point of Fig. [S] is that one sees quahtatively different trends for the scahng rela- 
tion depending on whether the multicomponent {—s^'^'^/k-Q or r) or effective one-component 
quantities {—si^/k-Q or t^s) are used. When the former are adopted, one finds that a single 
scaling curve describes all data for which e/k^T > 2 (i.e., the data for which rj follows 



the "dynamically normal" trend of increasing with e/k^T; see Fig. 2 of ref. l45l ). Only 
dynamically anomalous state points with e/k^T < 2 (i.e., for which rj decreases with in- 
creasing e/k^T) deviate from the scaling curve. This is analogous to the excess entropy 
scaling behavior observed in recent simulations of liquid water.— Specifically, deviations 
from excess entropy scaling were also observed for dynamically anomalous state points 
(where water's relaxation times decrease with increasing density). As can be seen from 
the insets of Fig. [5] (also in the data of ref. |45|), different trends emerge when the effective 
one-component quantities — sy/ZcB or t^s are adopted in the analysis, which obscure the 
aforementioned qualitative connection between the scaling behaviors of two different types 
of fluids that exhibit dynamic anomalies. 

IV. CONCLUSIONS 

This investigation highlights some significant problems that arise when characteriz- 
ing real-space, static pair correlations of moderately polydisperse colloidal fluids by a 
commonly-used effective one-component approach. Specifically, our analysis of a model 
fluid that mimics the experimental behavior of colloids with short-range attractions shows 
that trends based on the one-component treatment are qualitatively inconsistent with those 
of the full multicomponent description; e.g., the former spuriously indicates that pair corre- 
lations weaken with increasing attractions over a wide range of conditions. A key problem 
associated with the effective one-component analysis is an artiflcial, polydispersity-induced 
"softening" of the apparent static correlations-due to averaging over different particle 
sizes-which in turn causes the entropic consequences of structural ordering to be signifl- 
cantly underestimated. As an example of the implications of this issue, we demonstrate 
how excess entropy scalings for the shear viscosity of our model fluid show expected behav- 
ior only when the full multicomponent structural analysis is employed. Finally, we show 
that characterization of the partial radial distribution functions in concentrated suspen- 
sions of model colloidal fluid allows for an accurate estimate of the thermodynamic excess 
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entropy, a quantity that cannot be reliably estimated from an effective one-component 
analysis. 

What is the maximum level of polydispersity (e.g., variance in the particle diameter 
distribution) for which the aforementioned effective one-component description can accu- 
rately capture the structural properties of multicomponent fluids? Furthermore, how does 
the answer depend on the nature of the polydispersity? At present, the answers to these 
important questions are unknown. We are currently studying both issues, and we will 
report our findings in a future publication. 

On the theoretical and methods fronts, the present study demonstrates that both an- 
alytical results from statistical mechanics and Monte Carlo simulations can provide in- 
sights into some of the key polydispersity-related issues for the structural characterization 
of model colloidal fluids. Another important open question is how to carry out a simi- 
lar multicomponent analysis for data obtained via, e.g., optical microscopy experiments, 
where one is inevitably limited by inadequate sampling. We are currently exploring the 
prospects of using maximum likelihood estimation techniques together with liquid-state 
theory for accomplishing this, and we will report our findings in a subsequent publication. 
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Appendix: Exact results for low particle density 

In the limit of low particle density (p — ?■ 0), some implications of adopting the effective 
one-component description can be worked out analytically for a continuously polydisperse 
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SW fluid with pair interaction Vij{r 





oo r < (Tij 


Mr) = 


-e a,j < r < (1 




r>(l + A)(T 



\y. 



«j 



(20) 



and a Gaussian distribution of particle diameters with mean a and variance v'^a'^, i.e., 
p(c"i) = (27rt;^cr^)^^/^exp[— (o"j — o")^/2i;^o"^]. Note that, for this system, the probability 
distribution for the interaction diameter, /((Tjj) is a Gaussian with mean a and variance 
z/^cr^/2. In the low density limit, the PRDFs are known exactly, i.e., 







r < (Tij 



QtiV^) 



e^'"-"^ (Jij <r < (1 + A)o-ij 



1 r > (1 + X)aij 



(21) 



Substituting eq. [2T] into the continuous polydispersity limit of eq. HJ —s^'^'/k-Q = 
{p/2) J daijf{aij) J dr[gij{r) \ngij{r) - gij{r) + 1], noting that (afj) = a^{l + 3z/V2), and 
evaluating yields 



,(2) 



kspa^ 



2tt 



l + ^Ml+[(l + A)^-l] 



^e/ksT 



knT 



-1+1 



(22) 



The right-hand-side of eq. |22]is equal to B/a^ + [d{B / cx^) /d In T], where B is the second 
virial coefficient of the fluid. Thus, to leading order in density, we also have s^^' = s'^^, as 
expected. Likewise, substituting eq. [21] into the continuous polydispersity limit of eq. O 
r = y~^ J daijf{aij) f^" dy\gij{y) — 1|, and evaluating gives the following expression 



pi/3a 



-{l + A(e^/^«^-l)} 



(23) 



As can readily be seen, the right-hand sides of both eq. [22] and [23] are monotonically 
increasing functions of e/ksT. In other words, no structural anomalies appear in the 
dilute SW fluid according to the exact, multicomponent analysis. This is consistent with 
the recently reported low-density behavior^^ of a related single-component model. 

The effective one-component radial distribution ges{r) = J daijf{aij)gij{r) can also be 
obtained from knowledge f{<Jij) and gij{r) of eq. [JT] it is given by 



9es {rj 



1 + (1 - e^/^'B-) erf \ ^^"^ " ^\+ ^H + e^^'^^eri ^ ^/" " ' 



(24) 
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Substitution of eq. IH] into eq. E] and [7] allows determination of —s^q and t^s, respectively. 
Note that these quantities, unlike — s'-^^ and r, show a minimum at intermediate values 
of e/k-QT , giving the spurious impression that pair correlations anomalously weaken with 
attractions at low values of e/ksT. 

As a final point, we note that the single- component system with structure equivalent 
to that indicated by the effective one-component analysis will have the following pair 
potential, Vefr(r) = —kBT In ges{r). Polydispersity makes this potential considerably softer 
than the SW interaction between any two particles in the polydisperse fluid. 
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FIG. 1: Effect tliat the magnitude of the interparticle attraction ejk^T has on the static structure 
of the polydisperse SW fluid described in the text in the dilute (po"^ -^ 0) limit, (a) [Main panel] 
Comparison of the reduced two-body excess entropy, —s^'^'/pa^k-Q (identical to —s'^^/pa^kB in 
this limit), computed using eg. 1221 ^^ well its effective one-component counterpart, —s^^pa^kB, 
calculated using eq.[5]andl24i [Inset] Comparison of the reduced structural order metric, r/p^'^o", 
computed using eqJTTl as well its effective one-component counterpart, Tcs/p^'^a, calculated using 
eq. [7] and [Ml Effective one-component quantities spuriously indicate the presence of structural 
anomalies (attractions apparently weaken structure) for energies e/k^T less than the values 
indicated by the black dots, (b) [Main panel] The effective one-component radial distribution 
function gcs{i^), computed using eq. [Ml plotted versus reduced interparticle separation r/a. 
[Inset] The interparticle potential, Vesii") = —k-^T In gcs{r), of the effective one-component fluid. 
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FIG. 2: Partial radial distribution functions gsoj{r) describing correlations between the 30 and 
the j'*'' pseudo-components (j G [1, 60]) of the polydisperse SW fluid obtained via MC simulations. 
Results are shown as a function of reduced interparticle separation r/a and interaction diameter 
(o'so + '7j)/2. Color (red-to-blue) corresponds to magnitude of g30j{r) (high-to-low). Data are 
for a reduced density of po"^ = 1.05 and interparticle attractions of (top panel) e/ksT = 
(hard-sphere limit) and (bottom panel) e/k-^T = 2.5. 
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FIG. 3: Effect tfiat magnitude of tlie interparticle attraction e/k-^T has on tlie static pair cor- 
relations of the polydisperse SW fluid at a reduced density of pa^ = 1.05. Data obtained via 
MC simulations, (a)-(f) Partial radial distribution functions g30j{r) between 30*'* and the j^^ 
pseudo-components for j = 10, 20, 30, 40, and 50 (with corresponding characteristic interaction 
diameters a^oj = 0.90, 0.95, 1.0, 1.05 and 1.10, respectively) and various values of e/ksT. (g) 
Radial distribution function for the effective one-component description of the fluid, gcs{^), as a 
function of reduced particle center separation r/a for conditions used to generate the data for 
panels (a)-(f). 
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FIG. 4: Effect that magnitude of interparticle attraction e/k-QT has on different structural order 
metrics for the polydisperse SW fluid at a reduced density of pa^ = 1.05 Data obtained via 
MC simulations. [Main panel] Negative thermodynamic excess entropy, —s'^^/k-Q (diamonds), 
and its two-body approximation based on multicomponent, —s^'^'/k^ (squares), and effective 

(2) 

one-component, — s^g/fce (circles), descriptions of the system. [Inset] Translational structure 
metrics based on the effective one-component, Tes (circles), and multicomponent, r (squares), 
descriptions of the system. Curves are guides to the eye. 
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FIG. 5: Reduced shear viscosity rip~'^/^{mk-QT)^^''^ plotted versus various structural order met- 
rics based on static pair correlations for the polydisperse SW fluid: (a) [main panel] negative 
two-body excess entropy, —s'^'^'/k-Q and [inset] its effective one-component counterpart, —s^^/k-Q] 
(b) [main panel] translational structure metric, r, and [inset] its effective one-component counter- 
part, Tcfj. Structural data were obtained via the MC simulations of this study, and the dynamic 
data were extracted from the earlier molecular dynamics investigation of Krekelberg et al.— . 



Open symbols represent data for e/k-^T > 2 and filled symbols for e/k-^T < 2. 



conditions, shear viscosity anomalously increases with e/k^T (see Fig. 2 of ref. 



45|) 



^or the latter 



26 



